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OBJECTIVE — Type 2 diabetes and obesity are increasingly 
affecting human populations around the world. Our goal was to 
identify early molecular signatures predicting genetic risk to 
these metabolic diseases using two strains of mice that differ 
greatly in disease susceptibility. 

RESEARCH DESIGN AND METHODS— We integrated meta- 
bolic characterization, gene expression, protein-protein interac- 
tion networks, RT-PCR, and flow cytometry analyses of adipose, 
skeletal muscle, and liver tissue of diabetes-prone C57BL/6NTac 
(B6) mice and diabetes-resistant 129S6/SvEvTac (129) mice at 6 
weeks and 6 months of age. 

RESULTS — At 6 weeks of age, B6 mice were metabolically 
indistinguishable from 129 mice, however, adipose tissue showed 
a consistent gene expression signature that differentiated be- 
tween the strains. In particular, immune system gene networks 
and inflammatory biomarkers were upregulated in adipose tissue 
of B6 mice, despite a low normal fat mass. This was accompanied 
by increased T-cell and macrophage infiltration. The expression 
of the same networks and biomarkers, particularly those related 
to T-cells, further increased in adipose tissue of B6 mice, but only 
minimally in 129 mice, in response to weight gain promoted by 
age or high-fat diet, further exacerbating the differences between 
strains. 

CONCLUSIONS — Insulin resistance in mice with differential 
susceptibility to diabetes and metabolic syndrome is preceded by 
differences in the inflammatory response of adipose tissue. This 
phenomenon may serve as an early indicator of disease and 
contribute to disease susceptibility and progression. Diabetes 
59:2960-2971, 2010 




Type 2 diabetes and obesity are major causes of 
mortality and morbidity worldwide (1). Accord- 
ing to World Health Organization estimates, 
more than 1 billion adults are overweight and 
more than 200 million individuals have type 2 diabetes. 
The etiologies of obesity and diabetes are complex and 
created by interactions between environmental factors 
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(i.e., high caloric intake and reduced energy expenditure) 
(1,2) and genetic background (3-5). 

Inbred mouse strains are useful models for studying the 
role of the environment and genes in differential suscep- 
tibility to diabetes and obesity (6-14). In response to 
high-fat diet, aging, or genetic challenge, C57BL/6NTac 
(B6) mice become severely insulin resistant, hyperinsu- 
linemic, and diabetic, whereas 129S6/SvEvTac (129) mice 
are resistant to these challenges (7,8). Using intercross 
and F2 mice, our group and others have previously shown 
that the difference in disease susceptibility between the 
strains is inherited in a dominant fashion and linked to 
quantitative traits on at least three different chromosomes 
(8,15). 

In the present study, we compared gene expression 
profiles of B6 and 129 mice in different tissues, ages, and 
diets. For each comparison, we identified the differentially 
expressed network of genes between the strains using a 
novel variant of gene network enrichment analysis 
(GNEA) (16), an algorithm that integrates gene expression 
data together with protein-protein interaction networks. 
We subsequently intersected the results of the compari- 
sons between B6 and 129 mice to identify genes and 
pathways that were significant to each or all conditions, 
i.e., different tissues, ages, and diets. Of the numerous, 
differentially expressed subnetworks in the adipose tissue 
of B6 versus 129 mice; the most significant related to the 
immune system. Importantly, these differences were iden- 
tified even at 6 weeks, an age when the mouse strains were 
indistinguishable by fat mass or metabolic phenotyping. 
RT-PCR and flow cytometry confirmed higher expression 
of major inflammatory markers and higher infiltration of 
macrophages and T-cells in adipose tissue of 6-week-old 
B6 versus 129 mice. Weight gain associated with aging or 
high-fat diet further increased inflammation in adipose 
tissue of B6 but not in 129, mice. Taken together, these 
results demonstrate that measurable differences in the 
inflammatory milieu of the adipose tissue precede measur- 
able differences in insulin sensitivity in response to weight 
gain and contribute to the differences in diabetes risk 
between B6 and 129 mice. 

RESEARCH DESIGN AND METHODS 

C57BL/6NTac and 129S6/SvEvTac male mice were obtained from Taconic 
(Germantown, NY). Mice were maintained on a 12-h light-dark cycle with ad 
libitum access to tap water (reverse osmosis purified) and a chow diet 
containing 21% calories from fat, 23% from protein, and 55% from carbohy- 
drates (Mouse Diet 9F; PharmaServ, Framingham, MA). For the feeding 
studies, 6-week-old mice were submitted to low-fat diet (Rodent NIH-31M 
Auto: 14% calories from fat, 25% from protein, and 61% from carbohydrates 
[Taconic]) or high-fat diet (TD. 93075: 55% calories from fat, 21% calories from 
protein, and 24% calories from carbohydrates [Harlan Teklad, Madison, WI]) 
for 18 weeks before they were killed. In addition to fat content, these diets 
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TABLE 1 

Metabolic parameters in B6 vs. 129 mice at different ages and dietary conditions 

Chow diet (21% fat) Low-fat diet (14% fat) High-fat diet (55% fat) 

6 weeks old 6 months old 



B6 129 B6 129 B6 129 B6 129 



Body weight (g) 


21.9 ± 


0.6 


23.5 ± 


: 0.5 


45.5 


+ 


1.0 


30.4 ± 1.2 


41.5 ± 


1.4 


30.9 ± 


1.2 


46.4 ± 


1.1 


39.2 ± 


1.5 


Epididymal fat 




































mass (g) 


0.09 ± 


0.01 


0.21 ± 


0.02 


1.67 




0.14 


2.14 ± 


0.30 


1.62 ± 


0.17 


0.58 ± 


0.11 


3.39 ± 


0.33 


3.92 ± 


0.14 


Blood glucose- 




































random fed 




































(mg/dl) 


155.5 ± 3.1 


129.3 ± 2.9 


168.2 


+ 


12.3 


106.6 : 


± 3.5 


148.5 ± 7.4 


102.5 ± 3.4 


153.3 ± 6.7 


105.5 ± 3.7 


Insulin (ng/ml) 


2.26 ± 


0.48 


3.12 ± 


0.58 


92.14 


+ 


25.71 


3.96 ± 


2.04 


24.78 ± 


12.51 


1.29 ± 


0.32 


32.73 ± 


12.46 


3.58 ± 


0.41 


Leptin (ng/ml) 


7.14 ± 


0.70 


9.00 ± 


0.90 


76.18 


+ 


3.82 


16.67 : 


± 4.5 


72.50 ± 


8.06 


4.97 ± 


0.77 


82.38 ± 


7.41 


46.81 ± 


6.02 



Values are expressed as mean ± SE of 8 animals per group. 



differ in other components that could also influence the results of our analysis. 
For example, our chow diet contained fat from animal sources, whereas both 
low-fat and high-fat diets contained vegetable fat. These differences may have 
contributed to the fact that a few metabolic parameters displayed in Table 1, 
such as insulin levels, did not necessarily track with the percentage of fat in 
the diet. However, in this study and in others, these diets have been primarily 
used as independent paradigms to investigate the impact of the environment 
on weight gain of animal models, showing metabolic effects that closely 
resembled different characteristics of human obesity (8). Caloric-restricted 
animals were obtained from the National Institute on Aging. Caloric restric- 
tion was initiated at 14 weeks of age with a 10% decrease in calories, increased 
to 25% restriction at 15 weeks, and to 40% restriction at 16 weeks, which was 
maintained until the age of 6 months when the mice were killed. Protocols for 
animal use were reviewed and approved by the Animal Care Committee of the 
Joslin Diabetes Center and were in accordance with the National Institutes of 
Health guidelines. 

Gene expression analysis. RNA pooled from tissues of two or three mice 
was used for microarray analysis. Two to five Affymetrix murine chips 
U74Av.2 (Santa Clara, CA) were used per mouse strain, for each comparison 
of B6 versus 129 mice. Quantitative RT-PCR was performed using a SYBR 
Green-based ABI Prism 7900 Sequence Detection System (Applied Biosys- 
tems, Foster City, CA). For details, refer to the supplementary methods in the 
online appendix available at http://diabetes.diabetesjournals.org/cgi/ 
content/full/db 10-0367/DC 1. 

Gene network enrichment analysis. A novel version of Gene Network 
Enrichment Analysis (16) was applied to each comparison of B6 versus 129 
mice to identify gene networks (subnetworks of genes with edges connecting 
physically interacting protein products) with differential activation between 
these mice, and to test such networks for statistically significant, over- 
represented pathways and processes. Briefly, the algorithm was composed of 
the following steps (see supplementary methods for details). 

First, for each gene, a differential expression significance value was 
calculated in B6 versus 129 mice, converted to z scores, and used to annotate 
the corresponding proteins in a protein-protein interaction network curated 
from literature (17) (Fig. L4). 

Subsequently, a computationally intensive, stochastic optimization algo- 
rithm (18) was employed to identify a subnetwork of genes with a high mean 
z score (reflecting high significance). This process was repeated many times 
to yield a multitude of subnetworks (Fig. IE). The hypothesis was that a 
consensus among the stochastically detected subnetworks represented a 
prediction of the true, differentially active gene network in B6 versus 129 mice. 
Thus, each gene was associated with a P value based on the observed number 
of subnetworks in which it appeared, reflecting its chance probability of 
belonging to the true network (Fig. 1(7). 

In the next step of the algorithm, each biological process and pathway from 
the Gene Ontology and Molecular Signatures Database were tested for 
over-representation of genes with low P values, using standard statistical 
methods (Fig. ID). 

Finally, for each statistically significant pathway or biological process, the 
algorithm produced a gene network consisting of pathway member genes and 
all additional genes with significant, adjusted P values (false discovery rate 
<0.25) that shared an edge with a pathway member gene in the protein- 
protein interaction network (Fig. IE), providing an intuitively appealing 
visualization aid that helped in interpreting the enrichment results. 
Flow cytometry. Erythrocyte-free stromavascular cells were isolated from 
epididymal fat pads after collagenase digestion as previously described (19). 
These cells were incubated with a mix of antibodies against different surface 
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markers and analyzed using the LSRII flow cytometer (BD Biosciences, San 
Jose, CA). See supplementary METHODS for detailed protocol. 
Statistics. Statistical analyses for GNEA are described in supplementary 
methods. All other experiments were compared using ANOVA. Results are 
expressed as mean ± SEM. P values < 0.05 were considered significant. 

RESULTS 

Gene network signatures associated with the predis- 
position to metabolic diseases in B6 versus 129 mice. 

At 6 months of age, on a chow diet, B6 mice were 
overweight, hyperglycemic, and exhibited higher plasma 
leptin and insulin levels than 129 mice of the same age 
(Table 1, supplementary Fig. SI) (8). B6 mice showed 
clear signs of insulin resistance when compared with 129 
mice, even on low-fat diet (LFD), and this difference was 
further exacerbated when the mice were placed on a 
high-fat diet (Table 1) (8). At 6 weeks of age, however, no 
phenotypic differences or serum markers indicative of 
diabetes and obesity could be detected between the two 
strains (Table 1) (8). Indeed, 6-week-old 129 mice have the 
same weight as B6 mice and actually have significantly 
more fat mass, as indicated by higher fat pad weight (Table 
1). By studying B6 versus 129 mice at different ages and on 
different diets, one may therefore elucidate the back- 
ground inherited phenomena predisposing to obesity-ac- 
companied insulin resistance. 

To identify gene networks that might provide early 
biomarkers of adult metabolic disease, we initially per- 
formed microarray analyses on RNA from adipose tissue, 
liver, and skeletal muscle of B6 and 129 mice at 6 weeks 
and 6 months of age on a standard chow diet. We then 
applied GNEA to detect networks that were differentially 
active in a predisease stage and to estimate how they 
changed over time in disease-prone versus diabetes-resis- 
tant strains (Fig. 1). Since disease susceptibility is an 
inherited phenomenon in these animals, we hypothesized 
that the predisposing hereditable biological processes 
would be those showing differential activity in young mice 
and maintained or further increased with age. 

Among the three different tissues in which B6 versus 
129 mice were compared, only adipose tissue exhibited 
significantly enriched biological processes (Q value 
<0.33) at 6 weeks of age by GNEA, with 12 enriched 
gene sets, and this increased to 25 enriched gene sets at 
6 months of age. Of these, only two sets ("immune 
system process" and "regulation of cytokine secretion") 
were found enriched at both ages (supplementary Table 
1), and both of these related to inflammatory pathways. 
Furthermore, the "immune system process" was the 
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Experimental Setup 
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FIG. 1. Overview of gene network enrichment analysis. Individual comparisons of B6 and 129 mice were performed in different tissue types 
(adipose, skeletal muscle, and liver), ages (6 weeks and 6 months), and diets (high fat, low fat, and standard chow). A: For each comparison, a 
differential expression significance value was determined for each gene, converted to a z score, and mapped onto a network of protein-protein 
interactions. The colors in the network represent the different ranges of z scores. B: A stochastic search algorithm (18) was then applied to search 
for subnetworks with high mean z score (high significance). Each application of the algorithm identifies one subnetwork; hence, to detect a 
multitude of subnetworks, the algorithm was run many times. The consensus among these subnetworks represented approximations of the true, 
differentially active network in B6 vs. 129 mice. C: A membership in differentially active subnetwork (MAN) score was determined for each gene. 
MAN scores are P values and represent the probability of genes to belong to the true, differentially active network. This score was calculated as 
the random chance probability that the gene is detected in the observed number of subnetworks computed over permutations of the data. D: 
Biological processes and pathways were tested for enrichment in genes with high MAN scores. Statistically significant pathways and processes 
were predicted to be differentially active in B6 vs. 129 mice. E: For each significant process and pathway, a network visualization was generated 
to aid interpretation. 



most significantly enriched gene set at 6 months of age 
(Q value = 0.008; supplementary Table 1). In contrast, in 
the liver, 37 gene sets (Q value <0.33) were differen- 
tially regulated in B6 versus 129 at 6 months of age, but 
none was significantly different in the 6-week-old ani- 



mals. In skeletal muscle, no significant difference was 
observed at either age (data not shown). 

In addition to GNEA, the adipose tissue dataset was 
analyzed using gene set enrichment analysis (GSEA) (20) 
and the standard hypergeometric test (21). GSEA also 
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□ 6 weeks 
■ 6 months 




-4.85 -3.83 -2.82 -1.80 -0.79 0.23 1.24 2.26 
Fold-change (log 2 ) in B6 vs. 129 mice 



3.27 4.29 



5.31 



FIG. 2. Histogram of differential expression between B6 and 129 mice at 6 weeks and 6 months of age for genes in the immune system process 
network. The x-axis values are the upper bounds of the histogram bins and correspond to the log 2 fold-change difference in B6 vs. 129 mice. Few 
network genes are unchanged (absolute log 2 fold change <1) in B6 vs. 129 mice at 6 months. The majority are strongly upregulated (log 2 fold 
change ^1), although a few are strongly downregulated (log 2 fold change ^ -1). A network view of the same biological process demonstrating 
the fold changes of individual genes is provided in supplementary Fig. S2. 



identified differences in inflammation at 6 weeks. In par- 
ticular, it found a significant difference in the "humoral 
immune response" biological process (supplementary Ta- 
ble 2). On the other hand, GSEA failed to identify any gene 
set that was differentially expressed at both ages (supple- 
mentary Table 3). By comparison, the standard hypergeo- 
metric test (21) failed to identify either a clear immune 
signature at 6 weeks or any enriched gene sets at 6 months 
(supplementary Tables 2 and 3). 

Gene networks associated with the immune system 
process are differentially altered in the adipose tis- 
sue of B6 versus 129 mice. To determine how genes 
related to the immune response were differentially ex- 
pressed in the fat tissue of young B6 and 129 mice and 
changed in expression as the mice aged, we identified the 
network associated with the "immune system process" 
gene set (supplementary Fig. S2) and plotted the histo- 
gram of differential expression for its genes in B6 versus 
129 mice at 6 weeks and 6 months (Fig. 2). A few 
individual genes, such as apolipoprotein A (APOA1 and 
AP0A2) and mannose-binding lectin (MBL2) were more 
strongly downregulated in B6 versus 129 mice at 6 months 
and less so at 6 weeks. Nonetheless, much of the immune 
system process network was upregulated at 6 months 
compared with 6 weeks (note the right-shift in the histo- 
gram for 6 months vs. 6 weeks in Fig. 2). Lymphoid and 
myeloid markers, such as CD45, CD3, and chemokine (C-C 



motif) receptor 5 (CCR5), monocyte chemotactic protein-1 
(MCP1), as well as proinflammatory genes such as chemo- 
kine (C-C motif) ligand 5 (CCL5) and interleukin-6 (IL6), 
were more highly expressed in B6 compared with 129 mice 
at 6 weeks, and this increased further at 6 months of age 
(supplementary Fig. S2, left panel, arrows). Moreover, as 
shown by the histogram (Fig. 2), many of the immune 
system process genes already showed twofold or greater 
differences in expression in B6 versus 129 mice at 6 weeks, 
an age when no significant metabolic differences were 
observed between strains. 

Genes that have protein products that physically inter- 
act tend to correlate in expression (22) and are likely 
related by function or membership in the same signaling 
pathways. We consequently examined the subnetworks of 
genes that interact at the protein level with inflammatory 
markers identified in the strain comparisons (Fig. 3). For 
the T-cell marker Thyl, minor changes in the expression of 
its interacting genes were observed in the 6-week-old B6 
versus 129 mice comparisons (Fig. 3A, left panel). How- 
ever, by 6 months of age, more than 90% of the nodes in the 
network were differentially active between the two strains 
(large green triangles in Fig. 3A, right panel). Hence, genes 
interacting with Thyl were dramatically more upregulated 
with age-associated weight gain in adipose tissue of B6 
mice than in 129 mice, suggesting that B6 mice have more 
T-cells recruited to adipose tissue in response to increas- 
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FIG. 3. Gene networks associated with inflammatory markers in adipose tissue of B6 and 129 mice at 6 weeks (left) and 6 months (right) of age. Gene 
networks were generated by mapping genes that were significantly overrepresented (<? value < 0.25) among the gene network enrichment analysis 
results at 6 months intersected with those at 6 weeks onto protein-protein interaction networks involving at least one interactor to each inflammatory 
marker: (A) Thyl, (1?) CD45, and (C) MCP1. The CD45 network was further restricted only to genes showing a greater than twofold magnitude 
difference in B6 vs. 129 mice between 6 weeks and 6 months. This was done to aid visualization since the unrestricted network consisted of 416 
interactions and 141 genes. Colors range from bright red to green, corresponding to twofold less and twofold greater differences in expression between 
B6 and 129 at each age, respectively. Nodes point upward if the fold-change difference between B6 and 129 is greater at 6 months than at 6 weeks, and 
downward otherwise. The node size corresponds to the magnitude of that fold-change difference between the two ages. The inflammatory network 
around each biomarker is drawn to scale between the two ages. Networks around different biomarkers are not to scale with one another. 
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FIG. 4. Expression of inflammatory markers in 
adipose tissue of 6- week-old (6 weeks) or 
6-month-old (6 months) B6 and 129 mice. A: 
Leukocyte markers, including selective myeloid 
cell markers (e.g., macrophages and dendritic 
cells) and lymphoid markers (e.g., lymphocytes). 
B: Cytokines and chemokines implicated with 
inflammatory responses. C: Mediators of LPS 
signaling. mRNA expression was assessed by 
qPCR. All values are normalized by TATA-bind- 
ing protein (TBP) and expressed as fold change 
of the 6-week-old B6 average value. Results rep- 
resent mean ± SEM of 5-7 animals. *P < 0.05; 
**P < 0.01; ***P < 0.001 vs. 129. #P < 0.05; 
##P < 0.01; ###P < 0.001 vs. 6 weeks. 



ing adiposity than 129 mice. The same pattern was ob- 
served in subnetworks generated for CD45, a leukocyte 
marker (Fig. 35); CDllc, a dendritic cell marker (supple- 
mentary Fig. S3A); CD lib, a macrophage marker (supple- 
mentary Fig. S3B); MCP1, a monocyte chemokine (Fig. 
3(7); and tumor necrosis factor a (TNFa), a proinflamma- 
tory cytokine (supplementary Fig. S3C). Together these 
data suggest that B6 mice show a proinflammatory ten- 
dency in fat at 6 weeks of age not observed in 129 mice, 
and that this differentially increases with weight gain 
promoted by aging, affecting the pool of inflammatory 
cells. 

Changes in expression of inflammatory markers and 
cytokines/chemokines in adipose tissue of B6 versus 
129 mice. To validate the observations obtained using 
network analysis, we measured markers of inflammatory 
cells and cytokine/chemokine genes in the adipose tissue 
of B6 and 129 mice using quantitative RT-PCR (qPCR) 
(Fig. 4). In agreement with the microarray data, at 6 weeks 
of age, there were higher levels of macrophage/myeloid 
cell markers (CD68, CD lib, and CD 18), lymphocyte mark- 
ers (CD3, CD72, and CD80), T-cell attractants/activators 
CCL5/RANTES, stroma cell-derived factor-la (SDFla), 
and interferon-7 (IFN7), among others, in B6 versus 129 
mice (Fig. 4A and E). For most of these genes, the 
differences between strains were maintained or further 
exacerbated at 6 months of age. MCP1, a monocyte 
attractant chemokine, despite being similarly expressed in 
young B6 and 129 mice, exhibited a 4.4-fold increase in B6 
mice between 6 weeks and 6 months of age, whereas in 
129 mice, there was a 78% reduction with age (Fig. 4S). For 



a few cytokines, namely TNFa, IL6, and macrophage 
migration inhibitory factor (MIF), no differences were 
observed between the strains at any age (Fig. 4S). 

We also validated by RT-PCR the expression of selected 
genes identified as differentially present in the differen- 
tially active subnetworks identified by GNEA. Two medi- 
ators of lipopolysaccharide (LPS) signaling, i.e., the 
lymphocyte antigen 86 (Ly86) and LPS-binding protein 
(LBP), were highly ranked by GNEA and were substan- 
tially higher in B6 compared with 129 mice at both ages (Fig. 
4(7). 

To determine if the higher expression of inflammation 
markers in B6 versus 129 mice was restricted to the 
adipose tissue or present in other tissues, we assessed the 
expression of two macrophage markers (CD68 and F4/80), 
two T-cell markers (CD3 and Thyl), and two chemokines 
(MCP1 and SDFla) in the liver, skeletal muscle, and 
spleen of B6 and 129 mice at both ages using qPCR. The 
pattern of expression in liver was qualitatively similar to 
that in adipose tissue, but with smaller fold differences 
(supplementary Fig. S4A). In contrast, no differences were 
observed in the skeletal muscle or spleen in the expression 
for most of these genes (supplementary Fig. S4S and C). 
The two exceptions were MCP1, which was increased in 
all tissues in the 6-month-old B6 mice when compared 
with 129 mice, and CD3, which increased in both B6 and 
129 mice with age in skeletal muscle only (supplementary 
Fig. S4). 

An altered repertoire of inflammatory cells in the 
adipose tissue of B6 versus 129 mice. To determine if 
there was a cellular shift in the adipose tissue contributing 
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FIG. 5. Immune cell repertoire in adipose tissue of 6-week-old (6 weeks) or 6-month-old (6 months) B6 and 129 mice. Cells in the stromavascular 
fraction were labeled with fluorescence-conjugated antibodies against different myeloid cell and lymphoid markers and counted using flow 
cytometry. Nonerythrocytes stromavascular cells, PI-/Ter-119-; leukocytes, PI-/Ter-119-/CD45 + ; macrophages, PI-/Ter-119-/CD45+/CDllb+/F4/ 
80+/CD11C-; CDllc + -Myeloid cells, PI-/Ter-119-/CD45+/CDllb+/F4/80+/CDllc+; T-cells, PI-/Ter-119-/CD45+/CD3 + ; T-helper cells, PI-/Ter-119-/ 
CD45+/CD3+/CD4+.^4: Represents number of cells per gram of fat tissue and (2?) the total number of cells in the epididymal fat depot. All values 
are mean ± SEM of 4-8 animals. *P < 0.05; **P < 0.01; ***P < 0.001 vs. 129 mice. #P < 0.05; ##P < 0.01; ###P < 0.001 vs. 6 weeks. 



to our computational and qPCR results, we analyzed the 
inflammatory cell repertoire in the fat tissue of B6 versus 
129 mice using flow cytometry (Fig. 5). This revealed a 
significant twofold increase in the number of nonerythro- 
cyte stromavascular fraction cells per gram of epididymal 
fat in B6 compared with 129 mice at both ages (Fig. 5^4). At 
6 weeks of age, 75% of this difference was explained by a 
higher number of leukocytes in B6 fat versus 129 fat. 
Although the number of leukocytes per gram of adipose 
tissue decreased in both strains with age as adipose mass 
increased, at both time points the number of leukocytes 
was about 2.5-fold higher in the stromavascular fraction 
from B6 fat than in 129 fat (Fig. 5A). Likewise, the number 
of macrophages per gram of fat decreased with age 
because of growing fat mass, but was higher in the B6 
versus 129 mice at both ages (Fig. 5A). These data corre- 
lated with histologic evidence of crown-like structures, 
i.e., collections of macrophages surrounding dead or dying 
adipocytes, in the adipose tissue of 6-month-old B6 mice 
(supplementary Fig. SI). 

Somewhat surprisingly, the density of CDllc + -myeloid 
cells, which have been previously associated with insulin 
resistance states in mice (23), was decreased in B6 versus 
129 mice at 6 weeks and 6 months of age by 65 and 37%, 
respectively (Fig. 5A). In contrast, the T-cell population 
was dramatically increased by 7.1-fold in B6 mice at 6 
weeks of age (Fig. 5A). Furthermore, the density of T-cells 
increased by 25-fold in B6 mice with age-associated weight 
gain, but snowed only a threefold increase in 129 mice 
(Fig. 5A). Thus, at 6 months of age, B6 mice had > 55-fold 
more total T-cells and >60-fold more CD4 + T-helper cells 
in the adipose tissue compared with the 129 mice (Fig. 
5A). B-cells were also present in the adipose tissue of 
these mice; however, they were not considered in the 
analysis because of their modest numbers in comparison 
with other leukocytes (<1%). 

A similar pattern of lymphocyte/monocyte infiltration 
was observed in the 6-week-old mice when the total 
number of cells in the epididymal fat pads was calculated, 
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i.e., cells per gram X total grams of fat (Fig. bE). When 
expressed as total cells, however, it was apparent that the 
number of each of the myeloid and lymphoid cell popula- 
tion increased with age in both strains. Importantly, the 
differences between strains were maintained in both 
younger and older mice (Fig. 5B). These results confirm 
the gene network and expression analysis and demon- 
strate that most of these differences are consistent with 
differences in the lymphoid/myeloid population in the 
adipose tissue of B6 versus 129 mice. 
Effects of high-fat diet on the inflammatory status of 
the adipose tissue of B6 versus 129 mice. To investi- 
gate how differences in the genetic background could 
impact the effects of diet-induced obesity on the inflam- 
matory cell repertoire of the adipose tissue, cohorts of B6 
and 129 mice were placed on a high-fat diet (HFD) for 18 
weeks and compared with mice on LFD. On both LFD and 
HFD, the mRNA expression of the macrophage marker 
CD68 was higher by twofold or greater in B6 compared 
with 129 mice (Fig. 6A). Similarly, MCP1 expression was 
higher in B6 than in 129 on LFD, and this was maintained 
with HFD (Fig. 6A). The macrophage marker F4/80 and 
SDFla both trended higher in B6 than in 129 mice on LFD, 
but neither reached statistical significance (Fig. 6A). Inter- 
estingly, the HFD-induced increases in the expression of 
the T-cell markers CD3 and Thyl were completely abro- 
gated in the 129 mice (Fig. 6A). Caloric restriction, an 
intervention that drastically reduces fat mass, helps pre- 
vents age-associated diseases, and prolongs life span (24), 
decreased the expression of CD3 and Thyl in the adipose 
tissue of B6 mice by 41 and 51%, respectively, without 
affecting the expression of macrophage markers (supple- 
mentary Fig. S5). These results demonstrate a positive 
correlation between adiposity and the expression of T- 
cells in adipose tissue of B6 mice, a phenomenon that is 
significantly impaired in 129 mice. 

Accordingly, analysis of B6 versus 129 mice under LFD 
and HFD using network analysis revealed a substantial 
number of genes in the "immune system process" gene set 
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FIG. 6. Expression of inflammatory markers in adipose tissue of 6-month-old B6 and 129 mice fed with LFD or HFD. A: Expression of mRNA 
was assessed by qPCR. All values are normalized by TBP and expressed as fold change of the average value of the LFD B6 mice. Results 
represent mean ± SEM of 4-6 animals. *P < 0.05; **P < 0.01; ***P < 0.001 vs. 129 mice. #P < 0.05 vs. LFD. B: The network view of the 
immune system process differences between B6 and 129 mouse strains on a LFD (left) and HFD (right). The gene network was generated 
by mapping genes that were significantly overrepresented in B6 vs. 129 mice (Q value <0.25) among the GNEA results on HFD intersected 
with those on LFD onto protein-protein interaction networks involving genes annotated with the immune system process gene set. Genes 
in red are more than twofold higher in 129 mice compared with B6 mice; genes in green are more than twofold higher in B6 mice compared 
with 129 mice. Genes are denoted by a downward triangle if the fold change of the HFD animals is less than those on the LFD, and by an 
upward triangle for the converse; the size of the triangle denotes the magnitude of the difference. Arrows represent genes with specific 
interest commented on in the main text. 



with higher expression in B6 compared with 129 mice, as 
evidenced by the large, bright green nodes in Fig. 6B. The 
majority of these genes had greater differences between 
B6 and 129 mice under HFD than under LFD conditions, as 



shown by the high number of upward triangles in Fig. 6. 
In particular, T-cell marker CD3 and the subnetwork 
around T-helper marker CD4 were among the most up- 
regulated genes with HFD in B6 versus 129 mice (arrows 
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in Fig. 6B). Conversely, anti-inflammatory cytokines, such 
as interleukin-4 (IL-4) and interleukin-10 (IL-10) were 
reduced after HFD in the B6 mice, but not in 129 mice 
(arrows in Fig. 6B). 



DISCUSSION 

Although a relationship between obesity and inflammation 
in fat has been previously observed, in most cases these 
inflammatory changes have been viewed as being second- 
ary to obesity. Our data, comparing diabetes-prone B6 
mouse and diabetes-resistant 129 mouse at 6 weeks of age, 
show alterations in the inflammatory process in adipose 
tissue even before differences in metabolic parameters can 
be detected. Thus, B6 mice exhibit increased expression of 
the T-cell chemokines SDFla and CCL5/R ANTES and an 
increased number of T-cells in the fat tissue. These differ- 
ences are associated with higher IFN7 and CD80 levels in 
the B6 mice — both molecules are known to participate in 
T-cell function and activation (25) (summarized in the 
model shown in Fig. 7). 

Recent work has shown that T-cells infiltrate into the 
visceral adipose tissue of obese animals and humans with 
type 2 diabetes, and this is followed by recruitment of 
macrophages and development of insulin resistance (26- 
29). Based on these findings, a model for the role of T-cells 
in the pathogenesis of obesity and insulin resistance has 
been proposed in which increases in SDFla and CCL5/ 
RANTES levels in adipose tissue occur in response to an 
obesogenic environment and promote infiltration with T 
lymphocytes (29). IFN7 derived from these T-cells then 
promotes MCP1 secretion by preadipocytes (and possibly 
other cell types), resulting in recruitment of macrophages 
that further contribute to insulin resistance by production 
of proinflammatory cytokines (29). 

Our study provides evidence to support a major impact 
of the genetic background of different mouse strains in the 
migration of T-cells to the adipose tissue, both in the basal 
state and in response to weight gain. In B6 mice, the 
number of T-cells in the adipose tissue correlates posi- 
tively with the increase in adipose mass as a consequence 
of aging or HFD. By contrast, this response is practically 
absent in 129 mice despite a significant increase in adipos- 
ity in response to age or HFD. Thus, similar to mice with 
ablation of T-cells (27,28), 129 mice develop only mild 
insulin resistance in response to obesity. To what extent 
aging and the composition of the diet impact the infiltra- 
tion of T-cells into adipose tissue, in addition to the effects 
of weight gain, remain to be determined. It is clear though 
that insulin resistance in response to increased adiposity 
differ substantially between mouse strains, and this phe- 
nomenon is correlated with the migration of T-cell to the 
adipose tissue. Thus, inflammation in the adipose tissue 
does not always correlate with weight gain and is strongly 
dependent on the genetic background of the host. Similar 
differences dependent on genetic background are likely to 
occur in humans and contribute to differences in obesity- 
induced diabetes risk in different ethnic groups or even 
different individuals, allowing for some of the "fat-fit" 
phenotype. 

Among other implications, these observations clearly 
impact the choice of models for metabolic studies, in 
particular when knockout mice are used. Knockout mice 
are often derived from the embryonic stem cells of 129 
mice, and therefore, in many cases, studies are performed 
on mice on a pure 129 background or mixed 129/B6 
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background. In these cases, backcrossing to B6 mice for 
several generations usually potentiates the inflammatory 
and metabolic response of the model to environmental 
factors. However, it is not clear which of these back- 
grounds (B6, 129, or mixed) most closely mimic the human 
condition. 

Other recent reports support our findings that metabolic 
disease traits can be associated with alterations of inflam- 
matory gene expression networks in the adipose tissue 
and liver (30,31). However, to our knowledge, no study has 
previously demonstrated a measurable difference in in- 
flammation preceding any measurable phenotypic differ- 
ence associated with metabolic diseases. The novelty of 
our study is that, by comparing young B6 and 129 animals 
using a sensitive computational approach, we can focus on 
factors that can potentially predispose to disease in a 
prospective manner and avoid findings that are mainly 
secondary to obesity or metabolic differences. Our inves- 
tigation also examined differences between B6 and 129 
mice in multiple organs, allowing us to conclude that 
inflammation in adipose tissue, and to a lesser extent in 
liver, but not in skeletal muscle or spleen, is associated 
with the predisposition to insulin resistance. In addition, 
our computational analysis was able to set forth hypothe- 
ses that led to subsequent biologic validation experiments 
which provided further insight into the components of the 
immune system that may contribute to metabolic diseases 
(i.e., T-cell recruitment). 

In addition to inflammation, other phenomena have 
been shown to participate in the multifactorial pathophys- 
iology of insulin resistance. Alterations in insulin receptor 
levels and insulin signaling through IRS-proteins (32), 
induction of the unfolded protein response (33) and oxi- 
dative stress pathways (34), and changes in lipid (35) and 
amino acid metabolism (36) can all promote insulin resis- 
tance and contribute to the final phenotype. Many of these 
pathways act by producing post-translational modifica- 
tions of signaling proteins, such as phosphorylation or 
alterations in compartmentalization, which would not be 
detected as changes at the gene expression level. In this 
regard it is worth noting that in addition to differences in 
inflammation, GNEA analysis was able to identify other 
gene sets differentially expressed between the mouse 
strains in adipose tissue at 6 weeks of age (supplementary 
Table S2), including networks related to signal transduc- 
tion, protein secretion pathways, and glucose catabolism. 
Many of these pathways can interact with inflammatory 
pathways, and this crosstalk could represent an entry 
point to the manifestation of metabolic diseases. Thus, 
although changes in the immune response are definitely 
one of the factors that precedes and predicts the tendency 
of B6 mice to have greater insulin resistance than 129, it is 
unlikely that inflammation is the only predisposing risk 
factor associated with diabetes between these two mouse 
strains. 

In summary, it has been proposed that type 2 diabetes 
and obesity are diseases associated with an immune 
system that cannot cope appropriately with environmental 
threats (37-39). Based on this hypothesis, anti-inflamma- 
tory drugs, such as salicylates (40,41) and interleukin-1 
blockers (42), have been used to improve glycemia in 
individuals with type 2 diabetes. In this study, we demon- 
strate that pre-existing differences in the inflammatory 
milieu in metabolically active tissues may represent an 
important component of the genetic background as a risk 
factor to metabolic diseases. Thus, inflammation cannot 
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FIG. 7. Schematic model illustrating the potential causes and consequences related to the different repertoire of immune cells in adipose tissue 
of B6 and 129 mice. Solid arrows, secretion; dashed arrows, migration; dotted arrows, migration/differentiation. 

be viewed as only a mechanism that links susceptibility mation is an important early variable in the metabolic 
factors such as overfeeding, underactivity, aging, and response to environmental challenges and suggests sev- 
stress to metabolic diseases — it may also link these pa- eral potential targets for intervention, including LBP, Ly86, 
thologies to heritability. Our study indicates that inflam- SDFla, CCL5/RANTES, and MCP1. This provides new 
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strategies for reducing the epidemic of type 2 diabetes and 
metabolic diseases in spite of increasing obesity by attack- 
ing the variable genetic risk. 
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